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THE ZENITH ANGLE DEPENDENCE OF COSMIC RAY PROTONS 


CHAPTER I 
INTRODUCTION 

The composition and spectral distribution of the 
cosmic radiation is a function of atmospheric depth, mag¬ 
netic latitude and the direction of observation. Here one 
refers, of course, to the total radiation. A complete 
description of the radiation would give the intensity and 
spectral composition of each class of particle as a function 
of these parameters and the time. Such a complete descrip¬ 
tion is not possible with present experimental techniques. 

In successively refined experiments (1-7) conducted by this 
laboratory, a technique of proton-meson analysis of the radi¬ 
ation has been developed and the spectral composition of 
these radiations studied at sea level and at 3-* kilometers 
altitude. These studies were all for radiation of vertical 
Incidence at the point of observation. It then became 
desirable to undertake further observations in order to 
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find the change in composition and intensity with direction 
of observation at fixed atmospheric depth. Two programs of 
this kind have been carried out. In the work of Charbonnier 
(6), the spectral composition of the radiation of range 
greater than 15 cm. of lead was studied as a function of 
zenith angle. As was shown, the radiation defined in this 
way consists almost entirely of mu mesons. As a result of 
these measurements, it was possible to give a much more 
detailed experimental check of the theory of atmospheric 
development of the mu meson radiation than had been pre¬ 
viously available. The present work reports a similar study 
of the radiation of range less than 15 cm. of lead. This 
radiation contains, in addition to low energy mesons, a 
large number of protons. Precautions taken to exclude the 
soft component are described in Chapter II, Section 1. In 
the past, of course, many studies of the zenith angle depen¬ 
dence of the total hard and soft components have been made. 
The present work is the first study of Its kind for any part 
of the nucleonic radiation. As protons and mesons at lower 
altitudes are not thought to be closely related genetically, 
it was expected that quite different dependences on zenith 
angle would be found for the two types of radiation. As it 
turned out, the angular dependences of protons and mesons 
were found to be the same, within experimental error. As 
will be shown, this must be considered a fortuitous result 






3 


of the particular condition* of- observation. 

Charbonnier‘8 observations on mesons were easily 
interpreted because of the almost complete lack of nuclear 
interaction of these particles. Their distribution as a 
function of depth and angle can be interpreted entirely in 
terms of the well understood processes of decay and energy 
loss by ionization. Such is not the case for protons or 
other particles showing strong nuclear interaction. Cross- 
sections are such that on the average a dozen or more colli¬ 
sions will have taken place between the top of the atmosphere 
and the point of observation of 3.4 kilometers. These colli¬ 
sions involve pi meson production along with production of 
secondary nucleons at large angles to the direction of the 
primary. The development of this radiation involves a cas¬ 
cade process qualitatively similar to that of the soft com¬ 
ponent. This has been discussed by Messel, Janossy and 
others (8-20). However, almost complete lack of either 
experimental or theoretical knowledge of cross-sections for 
the many processes Involved causes the result to be one that 
can hardly be expected to be experimentally confirmed. In 
addition, a simplified one dimensional cascade has been 
considered so that the result cannot be compared with the 
present studies of angular dependence. It had been expected, 
because of the Isotropy of the higher energy proton radiation 
plus the large angle scattering in successive nuclear collisions. 
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that ihe observed proton radiation at lower altitudes also 
would oe almost isotropic., It was surprising to find that 
the intensity decreased by a factor of 3-2 In going from the 
vertical to a zenith angle of 45 degrees. This result in 
effect places an upper limit on the average deviation suf¬ 
fered in a nuclear collision. Rather than attempt a 
detailed cascade calculation based on questionable assump¬ 
tions regarding cross-sections, a simplified phenomenological 
calculation has been made to give a description of the observed 
proton zenith angle dependence. 

As a by-product of these observations, additional 
information on the proton spectrum has been obtained, as 
well as data on the meson zenith angle dependence for energies 
lower than those observed by Charbonnier. 



CHAPTER II 


THE EXPERIMENT 

1. Experimental Arrangement 

A sketch of the experimental arrangement is shown in 
Figure I. The entire assembly may be rotated about a North- 
South horizontal axis, perpendicular to the plane of the 
figure, so that the cosmic ray beam at any zenith angle up 
to 50° may be accepted. Anti-coincidence counters Al, which 
lie outside the acceptance beam defined by counters Cl and 
C2, serve to limit the number of cases in which shower events 
are recorded. Only events of the type Cl + C2 - Al - A2 - 
A3 are recorded. A one-centimeter thick absorber In the form 
of a half cylinder of lead placed above counter C2, plus the 
absorber represented by the walls of the cloud chamber and 
the coincidence counters establishes the minimum range and 
momentum of accepted particles. The 15 centimeters of absorber 
between C2 and A2 and 3» of course, determines the maximum 
range accepted. As discussed by Todd (5), such range defin¬ 
ition serves to resolve protons and mesons into distinct 
momentum bands, in the absence of scattering and nuclear inter¬ 
action. To the extent that these effects are ignored, then, 
the apparatus accepts mesons of momentum 0.09 to 0.34 Bev/c, 
and orotons in the momentum band 0.4 to 1.0 Bev/c. In other 
respects the apparatus is the same as that used by Garrison 
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(7). The Inclusion of the one cm. absorber over C2 has 
removed a small overlap between the maximum meson momentum 
accepted and the minimum momentum of proton acceptance 
present in his experiment. This modification is not so 
large as to prevent comparison of present with previous 
work. 

As shown in Figure I, anti-coincidence coverage 
below the 15 cm. of lead is quite complete. A double layer 
of counters has been used below the absorber to prevent 
particle leakage between the counters. Along the sides, a 
single layer was considered adequate, as in any case few 
particles will suffer such wide-angle scatterings. To 
prevent amplifier overloading from too high a counting 
rate, three anti-coincidence input amplifiers have been 
provided. One was fed from counters A1 while the tubes 
below the 15 cm. absorber were divided into two groups 
feeding the other two channels. The counter groups so 
defined could be separately scaled to detect defective 
counters. 

As in previous work, the 2.5 cm. of lead absorber 
above A1 serve to insure recognition of electron showers. 
Such showers either actuate one of the counters Al, in 
which case the event is not recorded, or are recognized by 
their appearance in the chamber. This will be discussed 




later. 




acceptance beam 



EXPERIMENTAL ARRANGEMENT 
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2. Experiments Performed 

Three sets of measurements were taken with this 
apparatus located at Climax* Colorado, elevation 3*3 km. and 
magnetic latitude 48° North and in continuous 24-hour opera¬ 
tion during the months of August, September, and October, 

1951* Initial conjectures as to the results of the zenith 
angle measurements led to the conclusion that appreciable 
differences in intensity would require rather large angles 
of observation, so for this reason, an angle of 45 was 
chosen. The three sets differed only in the zenith angle 
to which the apparatus was tilted, a measurement being made 
at vertical incidence, at 45° to the East and 45° to the 
West of the zenith. Since the magnet could be tilted in 
only one direction with respect to the trailer in which It 
was installed, the trailer was oriented to give an East or 
West zenith angle. In each case, a reference structure on 
the magnet was made level in order to obtain the same amount 
of tilt East or West. The direction of tilt was changed 
periodically so that all data in a given direction were not 
obtained at once. Time did not permit the taking of data 
at more than one angle to the East and West plus the vertical. 

For the vertical direction, a total of 3482 pictures 
was taken, of which 1387 were suitable for measurement accord¬ 
ing to the selection criteria to be described later. The 
West measurements provided 1126 measurable tracks out of a 
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total of 3962 pictures, the East 1304 out of 4488, for a 
total of 11,932 pictures. The magnetic field was reversed 
occasionally to indicate any asymmetry from this source. The 
data for vertical and East zenith angles are divided about 
equally between field direct and field reverse, but unfor¬ 
tunately the West data were all obtained with the field 
direct. 

3* The Equipment 

A brief description of the principal components of 
the equipment is Included for the sake of completeness only, 
as their detailed description has appeared before (4, 5 , 6 , 7 ). 
It should be mentioned that the use and constant improvement 
of this equipment by the above investigators during the three 
years preceding the date of this experiment was responsible 
for smooth, continuous, and dependable operation throughout 
the experiment. 

Figure I shows the essential features of the magnet 
assembly. The cloud chamber proper is centered between the 
magnet poles, a conical hole in the forward pole providing 
access for photography, and a hole in the rear pole providing 
for a tube connecting the sensitive volume of the chamber 
with the expansion diaphragm at the rear. Illumination of 
the ohamber was provided by four Sylvania type R4330 flash 
lamps, whose light was focused by cylindrical lenses into 
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the chamber, and by two 300 -watt Incandescent lamps for 
visual Inspection. These accessories, plus the necessary 
Geiger counters and absorbers, were bolted to the inner face 
of the rear magnet. The gap of 2i inches between the inner 
magnet faces, where the sensitive volume of the chamber, 
Geiger tubes, etc. were located, was maintained by steel 
spacers at the right and left sides of the magnet. These 
spacers also served to complete the magnetic circuit. 

The chamber was filled with enough 60-40 n-propyl 
alcohol-water mixture to saturate the volume, and finally 
filled with argon to a total pressure of 20 pounds per 
square inch above atmospheric. Expansion is accomplished 
by releasing compressed air from behind a neoprene diaphragm 
through two pop valves electrically connected in parallel. 

The two pop valves provide faster air release than the single 
pop valve used previously, while the electrically parallel 
connection insures that the two valves open simultaneously. 

The magnets provide a field of 8200 gauss with an 
energizing current of 800 amperes, and are cooled by oil 
flowing through the windings at a rate of 60 gallons per 
minute. The oil in turn is cooled in an oil-water heat 
exchanger. The cooling water was obtained from an artificial 
lake whose constant head and temperature provided ideal 
chamber operating conditions. Further temperature stabili¬ 
zation was provided by a valve controlling a bypass on the 
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heat exchanger. Temperature drifts were detected by the 
unbalance in an AC bridge, one leg of which was immersed in 
the oil system. The error signal thus derived was then 
applied, through suitable circuits, to the task of opening 
and closing this valve. This system assured temperature 
stability to 0.01° Centigrade. 

Appropriate electronic equipment provided for the 
following functions. 

a. Event Selection. Simple twofold coincidence of 
resolving time 4 microseconds is obtained using a Rossi type 
coincidence circuit. Its output is delayed by a 5 micro¬ 
second univibrator whose output feeds a second Rossi circuit 
in which anti-coincidence mixing occurs. An anti-coincidence 
pulse of ten microseconds duration then arrives at the mixer 
five microseconds ahead of the coincidence pulse, thus 
blocking its further progress. The anti-coincidence circuit 
has the desirable feature of being continuously retriggerable. 
Thus, anti-coincidence blanking persists for ten microseconds 
after the last pulse received, no matter how frequently they 
arrive. 

b. Programming. This unit contains power supplies 
including Oeiger voltage, a scale of 16 , with which each 
channel of the event selector, as well as the output of the 
anti-coincidence train, can be separately scaled, and mechan¬ 
ical registers for recording the C and C-A events. 
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c. Cloud Chamber Control. The desired event, selected 
by the event selector, and recorded by the program unit, then 
triggers this unit, which carries out the following sequence 
of operations: 

1. Expands the chamber 

2. Removes the clearing field voltage 

3. Records the expansion 

4. After a time delay of 0.06 seconds to allow time 

for droplet growth, fires the flash lamps 

5. Recompresses the chamber and reapplies the clearing 
field voltage 

6. Provides circuit sterilization for a period of 90 
seconds to allow time for the chamber to recover 

7. Move8 the film to the next frame 

8. Monitors all supply voltages 

9. Provides for manual operation of the chamber and 
camera rewind, and sterilization to avoid the 
recording of spurious data, e.g., in the event 
of failure of the magnetic field. 










CHAPTER III 


CURVATURE MEASUREMENTS AND THE DETERMINATION OP 
ABSOLUTE RATES 

The determination of momentum of a particle is accom¬ 
plished by measuring the radius of curvature of the track 
and applying the formula. 

-5 

p - 3Br x 10 

where p is momentum In Bev/c, B is the magnetic induction in 
gauss, and r is the radius of curvature in meters. With a 
magnetic field of 8200 gauss as used, p « r/4. Curvature 
measurements of nearly 4000 tracks required a compromise 
between inaccuracy and the consumption of time. The method 
used involves the comparison of the curvature of a track 
whose image is projected on the screen of a comparator, with 
that of a standard arc. The process whereby these arcs were 
prepared and a complete description of the measurement tech¬ 
nique have been presented in previous work (4,5). More 
recently, the optical system of the comparator underwent an 
overhaul, and the magnetic field reversals mentioned in 
Chapter II, section 2, were performed to disclose any curva¬ 
ture distortion Introduced by the overhaul. The data pre¬ 
sented here, consisting only of low momentum particles. 


:»*5SHSS£Si£':Sr" > iim 
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showed no such effect. However, the succeeding experiment 
by Charbonnier (6), In which high energy mesons were accepted, 
was far more sensitive to such an error and did indeed reveal 
the presence of one. Both the source and correction of this 
fault have been explained exhaustively by Charbonnier. The 
error amounts to introduction of an additional spurious curv¬ 
ature of 1/40 reciprocal meters in such a way as to increase 
the curvature of tracks appearing convex upwards in the com¬ 
parator. As described by the above authors, the preparation 
of the standard arcs precisely cancels out the error for 
tracks of this particular kind. For tracks of the opposite 
sign of curvature, however, the compensating error inherent 
in the standard arcs no longer compensates, but aatually 
accentuates the discrepancy to the extent of l/20 reciprocal 
meters, i.e., twice the distortion. Hence, the tracks appear¬ 
ing convex downwards in the comparator, and th*se include 
negative particles with-the magnetic field direct and positive 
particles with the field reversed, must have their measure¬ 
ments corrected by this amount. In Charbonnier's case, this 
correction was necessary. In this case, it is relatively 
unimportant, but it is felt that the correction should be 
made because of the systematic nature of the error. 

The apparatus used in this experiment does not record 
absolute rates. They can be Inferred, however, provided we 
make several assumptions. An accurate record was kept of the 



so-called C-A rate* which is the rate at which counters Cl 
and C2 provide simultaneous outputs (within 4 microseconds) 
unaccompanied by anti-coincidence counts. The cloud chamber 
does not respond to C A counts which occur during the 90 
second sterilization period provided to allow thermal recovery 
of the chamber after each expansion. Neither does each C-A 
count represent a single particle traversing the telescope, 
since shower protection is not complete. Identity of the 
nature of the event causing the count is learned on examining 
the photographs, which are placed in one of the following 
three categories. 

a. Single Tracks, Such a photograph is characterized 
by the appearance of a single track in such a direction as to 
indicate traversal of both C counters. The length of the 
track, as well as the occurrence of obviously pre- or post- 
expansion tracks Is Ignored. 

b. Showers. These are photographs showing two or more 
time-coincident tracks. 

c. Misses. Occasionally a photograph will show a com¬ 
pletely empty chamber. This Is also interpreted as a shower 
phenomenon. 

It is assumed that the C-A counts occurring during 
the 90-second sterilization time represent the same statis¬ 
tical distribution of events as those photographed. This 
means that the single track rate can be computed from: 
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Single Track Hate - (C-A Rat.) 

This single track rate then represents the total rate cor¬ 
rected for the 90-second recovery time and the occurrence of 
showers. The statistical uncertainty in this rate, as quoted 
in the data sheets, is the sum of the uncertainty in the C-A 
rate and the statistical uncertainty in the number of single 
tracks. However, not all single tracks are measured, since 
we have assumed that tracks whose length does not exceed 70# 
of the chamber diameter cannot be measured accurately. As 
compared with straight tracks, this criterion discriminates 
against curved tracks to the extent of about 3# at the lowest 
momentum measured (0.09 Bev/c). Above .25 Bev/c the discrim¬ 
ination is negligible. Thus we compute. 


Number of Single Tracks of Momentum p in dp per unit time 

(Number of Measured Tracks of mom. p in dp) fOA _ _ _i, 

- 1 - (Total" Numb e r of Measured Tlr acfe) - (Single Track Rate) 

If we divide this quantity by the product of the area of and 
solid angle subtended by the telescope, herein called the 
aperture, and the width of the momentum interval, dp, we 
arrive at the absolute rate in particles per Bev/c cm sterad 
hour. If we use the formula on page 13, p - r/4, then dp - 
dr/4. For convenience, a normalization factor can be defined 
for each set of data: 


Normalization Factor - 


(Single Track Rate)(4) 


(Total Number of Measured Tracks) x 
(Aperture) 
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Thus: 

(No. of Tracks of Mom. p in dp) 
Absolute intensity - * 

The only ill-defined quantity included above is the 

aperture of the telescope. Its uncertainty traces to the 

uncertainty in the sensitive volume of the Geiger tubes. 

From the size and separation of the counters, this coverage 

2 • 

was estimated to be about 1 cm steradian, a figure later 
corrected to 0.94 ( 6 ) by comparison with the very accurate 
determination of the vertical intensity made by Greisen as 
quoted by Rossi ( 25 ). 





CHAPTER IV 


DATA 

Analysis of the data from the three sets of measure¬ 
ments appears in Tables I - III. The direction of the mag¬ 
netic field was reversed occasionally for the purpose of 
detecting the presence of any systematic curvature distortion. 
As mentioned in Chapter III, such a distortion is present, 
necessitating separate treatment for the two kinds of data, 
field direct and field reverse. In the 45° West case, only 
one direction of field was used, so that here no separate 
treatment is necessary. 

Listed In the tables are the total number of frames 
inspected and the number of frames containing showers, misses 
and single tracks respectively. -The number of frames con¬ 
taining single tracks whose length exceeds the lower limit of 
70# of the chamber diameter is listed as ’’Total Measured." 

The latter are further broken down into positive, negative, 
and "Infinite, no sign," these latter being in the measurable 
class, i.e., long enough, but of such small curvature that no 
definite assignment of sign could be made. 

In one case (45° East, field direct), chamber operation 
was so poor for a short period that accurate measurements were 
impossible. However, the tracks photographed during this 
short time were Inspected and classified as to type of event 
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(single track, shower or miss, positive or negative). Thus, 
the figures for these classifications in Table II include 
these data, but the total measured entry and the tabulated 
momentum spectra do not. The inclusion increases the statis¬ 
tical accuracy of the quoted rates and hence (but only 
indirectly) the accuracy of the spectra. The normalization 
factor and the various rates have already been explained In 
the chapter on determination of rates. Finally, the tabulated 
momentum spectra appear. The various entries are connected 
by the following formula; 

Absolute Intensity - ( g°di°s hfcevlffi ( N °™- p * ctor > 

+ ( Abs. Intensity) 

~ (No. of Particles)^ 


- --afc - 1'ifc 23afe -e* 
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TABLE I VERTICAL CASE 


Direction of Observation 
Magnetic Latitude. 
Altitude 


Vertical 
48° North 
3=3 Kilometers 


FIELD DIRECT DATA 


Positives 

Negatives 

Infinites, No Sign 
Total Measured 
Total Single Tracks 
Running Time 
C-A Counts 


655 

252 

3 

910 

1142 

173=33 hrs 
4732 


Total Frames 
Average Yield of 
Single Tracks 
C-A Rate per Hour 
Single Track Rate 
Normalization 
Factor 


0.0495 


Positives 


Letter 

Group 


Number of 
Particles 


Radius 

Interval 

meters 

0.316 

O .345 

0.429 

0.470 

0.380 

0.480 

0.660 

0.980 

1.790 

4.150 


Average 

Momentum 

Bev/c 

0.048 

0.120 

0.213 

0.322 

0.425 

0.530 

0.665 

0.865 

1.200 

1.920 


Absolute Intensity- 
particles per Bev/c 
cnr sterad hour 

1.7 + 0.6 

8.2 + l.l 
10.3 + l.l 

9.7 + 1.0 

5.5 + 0.9 

7.1 + 0.9 

6.2 + 0.7 
4.1 + 0.5 

1.7 + 0.25 
0.6 + 0.09 


Negatives 


O .312 

O .330 

O .390 

0.420 

O .320 

0.390 

O .510 

O .710 

1.140 

2.150 

13.330 


0.048 

0.116 

0.200 

0.300 

0.390 

0.475 

0,585 

0.730 

0.955 

1.350 


1.7 + 0.6 
7.1 + 1.0 
11.2 + 1.1 
8.4 + 1.0 
2.0 + 0.6 
0 . 25 + 0.1 


1+1+ 
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TABLE I (continued) 


FIELD REVERSE DATA 


Positives 
Negatives 
Infinite, No Sign 
Total Measured 
Total Single Tracks 
Total Frames 
Running Time 
C-A Counts 


^26 

148 


477 

63? 

1404 

82.64 

2404 


hrs 


Average Yield of 
Single Tracks 
C-A Rate per Hour 
Single Track Rate 
Normalization 
Factor 


45 % 

29.1 + 0.6 

13.1 + 0.9 

0.110 


Positives 


Letter 

Number of 

Radius 

Interval 

Group 

Particles 

meters 

0-A 

9 

0.312 

A-E 

23 

0.330 

E-I 

44 

0.390 

X-L 

55 

0.420 

L-N 

29 

O.320 

N-P 

31 

O.390 

P-R 

31 

0.510 

R-T 

43 

0.710 

T-V 

31 

1.140 

V-X 

22 

2.150 

X- 

8 

13.330 


Average Absolute Intensity 
Momentum particles per Bev/c 
Bev/c cm sterad hour 


0.048 

0.116 

0.200 

0.300 

0.390 

0.475 

0.585 

0.730 

0.955 

1.350 


3.2 + 1.0 

7.7 + 1.6 

12.4 + 1.9 

14.4 + 1.9 
10.0 + 1.9 

8.8 + 1.6 

6.7 + 1.2 

6.7 + 1.0 
3.0 + 0.5 
1.1 + 0.3 


Negatives 


0-A 

A-E 

E-I 

I-L 

L-N 

N-P 

P-R 

R-T 

T-V 

V-X 

X- 


4 

O.316 

0.048 

27 

O.345 

0.120 


0.429 

0.213 

42 

0.470 

0.322 

12 

O.38O 

0.425 

5 

0.480 

0.530 

2 

0.660 

O.665 

2 

O.98O 

O.865 

4 

1.790 

1.200 

1 

6 

4.150 

1.920 


1.4 + 0.7 
8.6 + 1.7 
11.0 + 1.7 

9.8 + 1.5 
3-5 + 1.0 
1.2 + 0.6 
0.3 + 0.2 
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TABLE II EAST CASE 


Direction of Observation: 45 East of Vertical 
Magnetic Latitude: 48° North 

Altitude: 3-3 Kilometers 


FIELD DIRECT DATA 

Positives 4l6 Average Yield of 

Negatives 240 Single Tracks 36.5# 

Infinite, No Sign 3 C-A Rate per Hour 10.1 +0.2 

Total Measured 533 Single Track Rate 3.7 + 0.2 

Total Single Tracks 767 Normalization 

Total Frames 2099 Factor 0.0278 

Running Time 287.85 hrs 

C-A Counts 2894 


Positives 




Radius 

Average 



Letter 

Number of 

Interval 

Momentum 

Absolute 

Intensity 

Group 

Particles 

meters 

Bev/c 

particles per Bev/c 
cnr sterad hour 

0-A 

4 

O .316 

0.048 

0.35 

+ 0.2 

A-E 

25 

0.3*5 

0.120 

2.1 

+ 0.4 

E-I 

39 

0.429 

0.213 

2.5 

+ 0.4 

I-L 

51 

0.470 

0.322 

3-0 

+ 0.4 

L-N 

18 

0.380 

0.425 

1.3 

+ 0.3 

N-P 

37 

0.480 

0.530 

2.1 

+ 0.35 

P-R 

39 

0.660 

O .665 

1.6 

± 0.3 

R-T 

46 

0.980 

O .865 

1.3 

+ 0.2 

T-V 

34 

1.790 

1.200 

0.5 

+ 0.09 

V-X 

26 

4.150 

1.920 

0.2 

± °-°3 

X- 

14 

-- 

— 

-- 

-- 



Negatives 



0-A 

5 

0.312 

0.048 

0.45 

+ 0.2 

A-E 

35 

0.330 

0.116 

2.9 

+ 0.5 

E-I 

80 

0.390 

0.200 

5.7 

+ 0.7 

I-L 

59 

0.420 

0.300 

3.9 

+ 0.5 

L-N 

6 

0.320 

0.390 

0.5 

+ 0.2 

N-P 

2 

0.390 

0.475 

0.14 

+ 0.1 

P-R 

2 

0.510 

0.585 

0.1 

+ 0.07 

R-T 

- 

0.710 

0.730 

- 


T-V 

- 

1.140 

0.955 

- 


V-X 

- 

2.150 

1.350 

_ 


X- 

- 

13.330 

— 

- 
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Positives 
Negatives 
Infinite, No Sign 
Total Measured 
Total Single Tracks 
Total Frames 
Running Time 
C-A Counts 


TABLE II (continued) 
FIELD REVERSE DATA 


404 

238 

645 

795 

2378 

379.87 hrs 
3602 


Average Yield 
of Single Tracks 
C-A Rate per Hour 
Single Track Rate 
Normalization 
Factor 


33.4# 

9.48 + 0.15 
3-06 + 0.15 

0.0190 


Letter 

Group 

0-A 

A-E 

E-I 

I-L 

L-N 

N-P 

P-R 

R-T 

T-V 

V-X 

X- 


Positives 


Number of 
Particles 

5 

27 

70 

70 

48 

58 

56 

35 

18 

11 


Radius 

Interval 

meters 

0.312 

0.330 

O .390 

0.420 

0.320 

O .390 

0.510 

0.710 

1.140 

2.150 

13.330 


Average 

Momentum 

Bev/c 

0.048 

0.116 

0.200 

0.300 

0.390 

0.475 

0.585 

0.730 

O.955 

1.350 


Absolute Intensity 
particles per Bev/c 
cm sterad hour 


0.3 

1.6 

3.4 

3.2 

2.8 

2.8 

2.1 

0.9 

0.3 


± 0.2 
+ 0.3 
+ 0.4 
± 0.4 
+ 0.4 
± 0.4 
+ 0.3 
+ O.15 
± °-07 


Negatives 


0-A 

A-E 

E-I 

I-L 

L-N 

N-P 

P-R 

R-T 

T-V 

V-X 

3 

49 

71 

63 

26 

8 

7 
. 4 

2 

h 

0.316 

0.345 

0.429 

0.470 

O. 38 O 

0.480 

0.660 

0.980 

1.790 

0.048 

0.120 

0.213 

0.322 

0.425 

0.530 

O .665 

O .865 

1.200 

0.18 

2.7 

3.1 

2.5 

1.3 

0.3 

0.2 

0.08 

X- 

*f 

2 

4.150 

1.920 

— 


0.09 + 0.05 


+ 0.1 
± 0.4 
+ 0.4 
1 0.35 
+ 0.3 
+ 0.1 
± 0.08 
+ 0.04 


IJ 
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TABLE III WEST CASE 


Direction of Observation: 45° West of Vertical 
Magnetic Latitude; 48° North 

Altitude; 3.3 Kilometers 


(All 

data with 

field direct) 



Positives 

775 

Average Yield of 



Negatives 

336 

Single Tracks 

35% 


Infinite, N 0 Sign 

16 

C-A Rate per Hour 

10.9 

+ 0.2 

Total Measured 

1127 

Single Track Rate 

3-8 

+ 0.2 

Total Single Tracks 

1338 

Normalization 




Total Frames 
Running Time 
C-A Counts 


3835 
508.70 hrs 
5538 


Factor 


0.0135 


Positives 


Letter 

Number of 

Radius 

Interval 

Average 

Momentum 

Absolute Intensity 
particles per Bev/c 
cm 2 sterad hour 

Group 

Particles 

meters 

Bev/c 

0-A 

14 

O .316 

0.048 

0.6 + 

0.16 

A-E 

54 

0.345 

0.120 

2.1 + 

0.3 

E-I 

101 

0.429 

0.213 

3.2 + 

0.3 

I-L 

109 

0.470 

0.322 

3.1 + 

0.3 

L-N 

63 

O. 38 O 

0.425 

2.2 + 

0.3 

N-P 

57 

0.480 

0.530 

1.6 + 

0.2 

P-R 

90 

0.660 

O .665 

1.8 + 

0.2 

R-T 

96 

O. 98 O 

O .865 

1.3 + 

0.15 

T-V 

78 

1.790 

1.200 

0.6 + 

0.1 

V-X 

65 

4.150 

1.920 

0.2 + 

0.03 

X- 

31 

-- 

-- 

- 

- 


Negatives 


0-A 

12 

0.312 

0.048 

0.5 

+ 0.15 

A-E 

66 

0.330 

0.116 

2.7 

± O .35 

+ 0.40 

E-I 

108 

O. 39 O 

0.200 

3.7 

I-L 

96 

0.420 

0.300 

3.1 

+ 0.30 

L-N 

15 

0.320 

0.390 

0.6 

+ 0.20 

N-P 

8 

0.390 

0.475 

0.3 

+ 0.1 

P-R 

7 

0.510 

0.585 

0.2 

T 0.1 

R-T 

4 

0.710 

0.730 

_ 

_ 

T-V 

4 

1.140 

0.955 

_ 

_ 

V-X 

11 

2.150 

1.350 

_ 

_ 

X- 

9 

13.330 

— 

- 

- 



TABLE IV 


COMBINED VERTICAL DATA 1950-1951 


Direction of Observation; 
Magnetic Latitude: 
Altitude: 


Vertical 
48° North 
3-3 Kilometers 


Positives 2666 

Negatives 1052 

Infinite s No Sign 31 

Total Measured 3749 

Total Single Tracks 5126 
Total Frames 10521 

Running Time 833.02 hrs 

C-A Counts 19,878 


Average Yield of 
Single Tracks 
C-A Rate per Hour 
Single Track Rate 
Normalization 
Factor 


48 . 7 # 

24.8 + 0.17 
12.1 + O .25 


0.0129 


Positives 


Letter Number of 
Group Particles 


Radius Average 

Interval Momentum 
meters Bev/c 


Absolute Intensity 
particles per Bev/c 
cm 2 sterad hour 


O-A 52 

A-E I 99 

E-I 321 

I-L 419 

L-N I 89 

N-P 252 

P-R 309 

R-T 331 

T -V 235 

V-X 193 


0.316 0.048 
0.345 0.120 
0.429 0.213 
0.470 O .322 
0.380 0.425 
0.480 O .530 
0.660 O .665 
0.980 O .865 
1.790 1.200 
4.150 1.920 


2.1 + 0.3 

7.4 + 0.6 

9.7 ± 0.6 
11.5 + 0.6 

6.4 + 0.5 

6.8 + 0.5 
6.0 + 0.4 

4.4 + 0.3 
1.7 ± 0.1 
0.6 + 0.05 


Negatives 


0 -A 

49 

A-E 

189 

E-I 

332 

I-L 

286 

L-N 

62 

N-P 

25 

P-R 

17 

R-T 

18 

T-V 

16 

V-X 

21 

X- 

37 


0.316 0.048 
0.345 0.120 
0.429 0.213 
0.470 0.322 
0.380 0.425 
0.480 O .530 
0.660 O .665 
0.980 0.865 
1 • 790 1 .200 
4.150 1.920 


2.0 + 0.3 

7.1 + 0.6 
10.0 + 0.6 

7.9 + 0.5 

2.1 + 0.3 
0.7 ± 0.2 
0.3 + 0.09 
0.2 + 0.06 
0.1 + 0.03 
0.07+ 0.02 
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CHAPTER V 


DISCUSSION OF THE RESULTS 

For the vertical case, data of-a similar nature had 
been obtained previously In studying the form of the proton 
spectrum (7). Nevertheless, observations for this direction 
were repeated. In part because of alteration of the equipment 
(Chapter II, Section 1) and m part as a check on the previous 
work, and because It was felt best that for comparison between 
the various directions of observation, the data should be 
taken in close succession rather than a year apart.. The 
present and previous spectra for the vertical case agree well 
within statistics, a fact that is regarded as proof that the 
equipment had undergone no change in performance during the 
intervening year . For improved statistics, the two sets of 
data have been combined with the result shown in Figure VI. 

An explanation of the difference in form of the spectra 
for positive and negative particles has been given previously 
in connection with the determination of the proton spectrum 
(7). It should be mentioned here that the negative distribu¬ 
tion contains mu mesons only, while the positive distribution 
is for positive mu mesons plus protons. While the total 
intensity of positive mesons is known to exceed that of the 
negative by some 20% s at momenta above 0.5 Bev/c where the 
selection method excludes almost all mesons the difference 
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between the positive and negative intensities can be taken as 
a measure of the proton Intensity. 

1. The Meson Spectrum 

. o 

In Figures II to V the vertical, the 45 East, and the 
45° West spectra are given, along with the 45° spectrum aver¬ 
aged over East and West The latter spectrum. Figure V, Is 
represented by plotted points only, the curves appearing on 
the figure for comparison purposes being the combined vertical 
data of Figure VI. The plotted points are the measured rates 
multiplied by a scale factor of 3.2 to allow comparison of the 
spectral shapes. With this multiplication factor and within 
statistics, the average East-West spectrum and the vertical 
spectrum do not differ. This is true for both the proton and 
meson radiation. An explanation of the zenith angle dependence 
of the meson radiation has been given in Charbonnler's study 
of this component. It was surprising to find the same zenith 
angle dependence for protons as for mesons since the two types 
of particles are not thought to be closely related genetically 
at this altitude. It should also be pointed out that the East 
and West spectra agree between themselves within statistics. 
This is shown in Figure VII. Thus, within the accuracy of 
the present observations, neither protons nor mesons at this 
altitude show any East-West asymmetry. The observed agree¬ 
ment In the zenith angle dependence of protons and mesons Is 
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FIGURE IV 
MOMENTUM SPECTRA 

DIRECTION-45* WEST 
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FIGURE V 
COMPARISON OF 
MOMENTUM SPECTRA 
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FIGURE VI 

MOMENTUM SPECTRA 
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FIGURE VII 
MOMENTUM SPECTRA 



MOMENTUM BEV/C 



34 




considered to be fortuitous. Thus while the meson intensity 
decreases by a factor of 32 from zero to 45 degrees for the 
low momentum mesons observed here, this ratio does not hold 
at higher energies. At 2 Bev/c, the intensity ratio for 
mesons is found to be about 2 (6). 

It is not the primary purpose of this experiment to 
explain the behavior of mesons. The zenith angle dependence 
of these particles has been thoroughly discussed by Charbon- 
nler (7). One assumes that essentially all mesons observed 
at lower altitudes are produced isotropically at great atmos¬ 
pheric heights. Taking Into account the negligible scattering 
of these particles, their radioactive decay and ionization 
energy loss In traversing the atmosphere, one accounts quan¬ 
titatively for the observed Intensity and directional distri¬ 
bution at lower altitudes. With these assumptions, a 
transformation can be found relating the vertical and 45° 
spectra at lower altitudes. In this way, from the observed 
vertical intensity at 3«3 km. for a momentum of 0.8 Bev/c 
one predicts a 45° intensity of 6.0 + .3 particles/Bev/c cm 2 
sterad hour, at 0.25 Bev/c. The 45° Intensity at 0.25 Bev/c 
observed In the ©resent work Is 7.19+0.33. Agreement 
between the predicted and observed values is not considered 
too bad, although outside the limits of statistical accuracy. 
Charbonnier did not observe meson intensities at momenta 
below 0.5 Bev/c. The predictions of a 45° intensity at 0.25 
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Bev/c was nevertheless made since this depends on his observed 
Intensity at 0.8 3ev/c. At this low residual momentum, because 
of the increased importance of scattering, the computation is 
not expected to have comparable validity with Charbonnier's 
results at higher momenta. As the radiation intensity 
decreases with increasing zenith angle, the net effect of 
scattering will be to increase the intensity at large zenith 
angles at the expense of the Intensity at smaller angles. 

This is in the direction of the observed and predicted 
intensities. 

2. The Proton Spectrum 

As compared to the meson radiation, the atmospheric 
development of the proton radiation (and the nucleonic radi¬ 
ation in general) is greatly complicated by nuclear interac¬ 
tion which here becomes the moat important single process. 
Development of the nucleonic radiation involves a cascade 
process qualitatively similar to that of the soft component. 

The cascade is propagated through the atmosphere principally 
by protons and neutrons (22-24) produced in successive pene¬ 
trating showers. Pi mesons originating in these showers 
decay too quickly to have high probability of undergoing 
further nuclear interactions while the mu mesons secondary 
to them do not undergo appreciable nuclear interaction at 
all, and so take no part in further propagation of the 
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cascade. Compared to the electromagnetic cascade, the inter- 

2 

action path lengths here are in the order of 100 g/cm , 

p 

instead of the few g/cm characterizing the radiation lengths 
of the electron cascade. The nucleonic cascade had been 
discussed in detail in a number of papers by Heltler, Janossy 
and Messel (8-20). The treatment there is phenomenological 
in regard to the form assumed for the nuclear cross-section. 
With this exception the mathematical treatment is complete 
with respect to fluctuations resulting from geometrical and 
statistical factors. The cross-section is assumed to be a 
function of the fractional energy loss of the primary, inde¬ 
pendent of the initial energy. The cross-sections employed 
do not include a description of the angular distribution of 
the secondary nucleons. Thus the zenith angle effects have 
not been computed, so that this treatment cannot be compared 
with the present experiment. 

The constants in the cross-sections used by the above 
authors were so chosen as to conform to experimental values 
for the removal path length of the nucleonic radiation. 
Absorption of the radiation producing penetrating showers 
has been studied by Van Allen, Bernardini, Rossi and others 
(25-29,33) who find an exponential absorption with a coeffi¬ 
cient of about 120 g/cm^ of air. Studies of the absorption 
of the radiation producing photographic stars as well as 
direct measurements by this laboratory of the relative proton 






intensity between 3.3 km and 3ea level gives similar absorp¬ 
tion path lengths. It should be pointed out that an expon¬ 
ential absorption of this path length cannot hold to the top 
of the atmosphere, as in this case the observed proton 
intensity of 3-3 km, would be much too large as compared to 
the primary proton intensity at the top of the atmosphere. 
Instead (and there Is some experimental evidence to the 
effect) it is probable that for the very high energy protons 
at the top of the atmosphere the effective removal path 
length Is much larger due to the larger number of "genera¬ 
tions" of secondaries of sufficient energy to continue the 
cascade. In the following treatment of the angular spread 
of the proton radiation a constant removal path length of 
120 g/cm is assumed. In making comparisons with observa¬ 
tions made at 3-3 km., it is not necessary that this 
assumption hold to the top of the atmosphere, but rather 
that it be good from heights of 3-3 km. to sea level. 
Protons observed in this Investigation lie principally 
between 0.5 Bev/c and 1.0 Bev/c. On the basis of ionization 
loss alone, these protons must have had momenta greater by 
some 2 Bev/c at the top of the atmosphere. As they will 
have undergone several nuclear collisions, the primary 
energy must have been several times greater than this. For 
protons of momenta above 7.5 Bev/c the primary distribution 
is isotropic above the horizon, as shown by the Lemaitre- 
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Vallarta theory (35) of the geomagnetic effect. It can be 
expected then that tne protons observed at 3°3 km. are 
derived from an initially isotropic radiation. The radiation 
at 3-3 km. then can be considered as the superposition of the 
radiation due to equal contributions from each direction of 
primary incidence. The observed angular spread at 3-3 km. 
as well as absorption path lengths as measured for either the 
total or vertical radiation depend in part on geometrical 
factors resulting from this superposition of contributions 
from many primary directions. It is the object of the fol¬ 
lowing treatment to attempt to separate out the purely 
geometrical factors and so arrive at an estimate of the 
angular spread and absorption pat-. length as determined by 
the physical processes and cross-sections. It is clear, for 
example, that a primary radiation all of vertical Incidence 
would not ^ive rise to so great an angular spread at 3=3 km. 
as that observed. Thus the observed angular spread repre¬ 
sents an upper limit to that due to nuclear scattering in 
traversing this absorber thickness. We attempt to obtain a 
closer approximation to the angular spread due solely to the 
nuclear scattering through the following assumptions: 

1. The primary proton Intensity is the same at all 
directions of Incidence. 

2. The intensity of secondaries at depth X due to 
primaries of zenith angle (independent of 



azimuth) depends only on the diagonal distance 

p 

d In g/cm as measured along this zenith direction. 

3. This intensity decreases exponentially with a path 

. 2 

length of 120 g/cm . 

4 . The angular spread of secondaries due to the primary 
radiation in a given direction is represented by a 
cos n id law centered about this direction with n 
independent of depth. 

The first assumption is completely justified for primary 
momenta above 7.5 Bev/c by the geomagnetic theory. As a 
result of the following computation the assumption that the 
primary radiation is isotropic for angles greater than 45 ° 
is found unnecessary. For this smaller angle, geomagnetic 
theory insures isotropy for momenta as low as 3-3 Bev/c. 

The second assumption would be strictly true if one 
neglected scattering. With scattering, the geometrical situ¬ 
ation is not strictly the same for all zenith angles. A 
primary traversing a given path near the top of the atmosphere 
will produce secondaries at lower altitudes with large lateral 
displacements from its original path. As compared to vertical 
incidence, minimum path lengths for these secondaries (from 
the point of observation to the top of the atmosphere) will 
be increased for 180° of azimuth and correspondingly decreased 
for the other 180°. Except for large zenith angles, since 
average path lengths are the same, the situation will not 
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differ from that for vertical incidence. 

Assumption three conforms to the large body of experi¬ 
mental determinations of this attenuation. While, as has 
been pointed out, this cannot hold to the top of the atmos¬ 
phere, for the present purpose it need only be a close 

approximation for atmospheric depths from about 700 to 1000 
. 2 

g/cm . Experimental evidence indicates its validity below 
60,000 feet (32). 

In regard to assumption four, a cos 1 ^ scattering law 
has no theoretical foundation. Its use conforms to usual 
practice, n being adjusted to give the experimentally deter¬ 
mined rate of decrease. The assumption that n is constant 
cannot be true over large ranges <f distance. Again this 
only need be approximately true over a range of 700 to 1000 
g/cm . As n will be a function of the proton energies 

involved, similarity of the proton spectra at 700 and 1000 
2 

g/cm can be taken as a measure of the constancy of n. 
Evidence for this has been provided by Wilson (21). 

The physical problem Involved here is similar to that 
considered In the Gross transformation (31)# with the addi¬ 
tion of the effects due to scattering. The additional 
complexity of the present situation necessitates resort to 
graphical methods. 

Symbols to be used are defined as follows (see Figure 
VII for the coordinate system used) . 
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I 0 {&o9^ 0 ) Differential proton Intensity at the top 

of the atmosphere at spherical coordinates 

2 

6 a m partlcles/Bev/c sterad cm 

hour.. This Is assumed to be constant in 
6 o and Vo > 

Io Total differential intensity at the top 

of the atmosphere <* /X 0 (9 0 <&) d-Oo « 

Jc^ ’ 

Z'n X a (&,<&) particles/cm 2 /Bev/c hour. 


x 


Depth of the point of observation measured 

2 

in path lengths of 120 g/cm . 


i A6o>%>9* , <f x ) 


i* 


Portion of the differential proton 

Intensity at depth x observed at angles 

O k s contributed by the primary 

radiation incident at angles 0 o » 0 at 

the top of the atmosphere, in particles/ 
2 2 

Bev/c cm sterad hour. 

Differential proton intensity at depth x 
observed at angles 0 X , due to the 


contributions from all primary radiation, 
in particles/Bev/c cm sterad hour. 


CO 


Angle between cUl© andd-Q*. See Figure 


VIII. 


i 


FIGURE VIII 

SCHEMATIC DIAGRAM OF THE 


COORDINATE SYSTEM USED 
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The counting rate of a Geiger counter telescope at 
the top of the atmosphere will be proportional to I (&, <f. )dOo 
where da 0 is the aperture of che telescope. The telescope at 
depth x will record I x (# x > We wish to account for 

the fact that this latter rate falls by a factor of 3.2 from 

o 

the vertical to 45 ' zenith angle. For this purpose we assume 
that ^ 

ija,,s’,) = kio(.9o,f.) Cos *oj 

We quickly find from Figure VIII that 

Cos co - 'S'n &o Cos (^ 0 - ^x) ■+• Cos 0 O Cos <0* 

so that 

= f kC 0 (& 0 , Vo) e c °* Go Q o Gxse K +5™o o Ci"0K Cos (<?«,- j] 

±o 

where k is a normalizing factor and n Q is the upper hemis¬ 
phere, The determination of k can be made immediately by 

noticing that f X* } <{ x ) J.Q x 

_>L_ r „ 

•= klo (Go, ?<>) e j Cos cj clSl K dSl 0 

J Slx 

must reduce to I Q % )dI2« for the case of no absorption, 
l.e.. Infinite removal path length. This means that x Is 
zero and the exponent!- 1 factor Is unity. Thus from 

k f [Cos Go Cos -I- S/» fifCos - ‘fx )1 d&x * I 

A* 
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where J 2 * is again the upper hemisphere, k can presumably be 
found. The result, however, evidently depends on 0 O • This 
is to be expected, since for larger values of G 0 , the cos n co 
distribution places particles outside the limits of integra¬ 
tion. However, it will be shown immediately below that the 
exponent n is large enough so that the contributions to the 
above integral for values of Ok greater than 90° are negli¬ 
gibly small for Go less than 45°, which is the region of 
interest. Thus it is legitimate to set 0 Q ■ 0, with the 
result * „ 

*/ / Cos”e^Sm0^<e K de x = / 

o o 

r f 

Ces n 0 x d(Cos0 x ) - / 


or 


k - 


A7 -h / 
Z 7T 


The determination of n is independent of this result. 

Involved in the choice of the cos n a> dependence for an 
assumed scattering distribution is the fact that forward scat¬ 
tering predominates and that scattering in excess of 90° is 
negligibly small. Since cos n <o vanishes at 90°, the behavior 
of the forward scattering is fairly well approximated. How¬ 
ever, cos n w is not zero beyond 90°, and contributions to the 
integral will occur unless appropriate limits for the inte¬ 
gration are defined. Hence, the integration will extend only 
ovr that port ion of the uppar hemisphere containing 
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directional rays making an angle less than or equal to 90 U 
with the direction defined by , Tt will be shown later 
that this limitation 13 unimportant. The physical situation 
demands that we have axial symmetry about the vertical (no 
deflection by the earth’s magnetic field is assumed ) , and 
hence, the result must be independent of This allows 

the choice - 0 - Thus the region of integration depends 
only on 6 ^ y and is shown in Figure IX. The shape of the 
region of integration shows that if the coordinates @0 , 
are retained, the limits of integration over ^ will be 
functions of <9 To avoid this complication, a change of 
coordinates is m order. This is shown in Figure X. The 
transformed Integral is then of ne form. T x (fi.<e x )dn K 

~ ,<&)[ Cos*(\s, n 5 Sin n +\ c/xdei dl 2 * 

'Jo 

where o( = — — Q» 

2 a 

We have data on this quantity for two directions of 
observation, hence the integration will be performed for 
just these two values, as it must be carried out graphically. 

The object of the computation, of course, is to 
find n. The values of I x (0 ><?* ) and I x (f><f x ) have been 
obtained for n « 5, 6 and 7 at x *= 5-4 corresponding to an 
atmospheric depth of 650 g/cm^ » 3-3 km. altitude. These 
are shown in Table V together with the ratios •* I J \O i ^ ) 















/l x ( ^- # ^ ) which are the measured quantities under discus¬ 
sion., It is seen that ^ Is a fairly sensitive function of 
n p and that by inrerrdating, r, can be assigned the value 6.5. 
Also included In the table is the value of for n equal to 
” infinity/’ i.e., assuring no scattering whatsoever and that 
is simply given by ^ where x' is the dis¬ 

tance to the top of the atmosphere along the diagonal ray at 
4 5 0 ;’,erlth angle to the vertical. 

The value of n so determined Is dependent on the choice 
of path length. The sensitivity of the dependence is not 

great. The figure for the path length as determined by this lab- 

/■* 

oratory Is 125 ± 8 g /cm;*. Other values (29) are more closely 

2 

centered about 120 g/cm . Values from nuclear plate data (26) 

2 

are somewhat, higher (about 135 g/cm ), but refer to particles 

in much broader momentum bands. Within the region L * 100 to 
2 

140 g/cm , n varies from 6.0 to 8.5. Within the region L » 

Q 

125 + 8 g/cm , and within the statistical accuracy of the 
experimental value of , n has the value 6,5 + 0„7. 

The value - 3.2 is thus shown to be the result of 
an effective scattering distribution of cos^‘(L . This is a 
fairly sharp distribution, falling off by a factor of nine 
from zero to 45 degrees zenith angle. It is of no small 
interest then that the use of such a distribution has the 
effect of reducing the value of p by a factor greater than 
three. The effect of scattering, however, must be included, 
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since the observed zenith angle dependence could be obtained 
purely from attenuation considerations only for path lengths 
in excess of 300 g/cm 2 . 

The behavior of the integrand of the above integral 
shows that the contribution to the vertical radiation at 
depth x from the primary radiation incident at angles greater 
than 45 degrees is entirely negligible (less than 0.1# for 
all values of n greater than 5 ), and contributions from the 
same source to the radiation observed at 45° at depth x is 
quite small (less than 5 # for n greater than 5), primarily 
as a result of the exponential factor. This provides the 
justification for the assumed isotropy of the primary radia¬ 
tion, since lack of isotropy at zenith angles greater than 
45° will have no effect on the computed values of I x (0 , ( f x ) 
and I x (^i <f x ) • Hence, we need justify an assumed isotropy 
only at angles less than 45° within which protons are iso¬ 
tropic at momenta as low as 3-3 Bev/c, according to the 
Lemaltre-Vallarta theory. 

Knowledge of the normalization factor allows one to 
compute an effective path length for removal which takes into 
account the scattering factor. Thus the computed ratio 


1 0 (o. <e») 

Wo. f*> 

(both in the vertical direction) represents the actual measured 
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p 

ratio of intensities at 650 and zero g/cm depths respec¬ 
tively, and when equated to exp(-650/L) yields the measured 
value of L. These ratios are exhibited in Table V for n *«- 5, 
6,7 and infinity, as well as the associated values for L. 

These effective values of L are all somewhat smaller than 
2 

120 g/cm , as is expected, since more particles are scattered 
out of the vertical beam than into it, thus increasing the 
apparent attenuation of the vertical beam. 


TABLE V 

INTENSITY RATIOS AND EFFECTIVE PATH LENGTHS 


n 


( 0 

,v 


*x(?.<e 

_ lx(0, <ftc) 

"Effective" 
Path Length 

5 

2.28 

X 

10~ 3 

*0 

.857 X 10 -3 I 0 

2.66 

IO 6.5 gm/cm' 

6 

2.46 

X 

io“ 3 

I© 

.840 x 10~ 3 I 0 

2.93 

108 

7 

2.65 

X 

io‘ 3 

I 0 

.755 x 10' 3 I o 

3.31 

109 " 

Inf. 

4.53 

X 

10' 3 

I© 

.465 x 10~ 3 I o 

9.70 

120 " 


As mentioned before, the exponential factor is respon¬ 
sible for the sharp cutoff at angles in the vicinity of 30-45 
degrees. As one goes to greater and greater depths, i.e., to 
larger values of x, the sharpness of the cutoff will be 
Increased somewhat (in addition to the larger attenuation) 
so that at sea level the result of a computation of the ratio 
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of the vertical intensity to the 45° intensity is about 3 . 5 . 
There is no experimental verification of this surmise,. The 
extreme difficulty of picking up any protons at all at sea 
level (about \% of the total penetrating component) indicates 
that a statistically meaningful determination of this ratio 
would require a long and careful but not impractical measure¬ 
ment . 

As mentioned in the introduction, the above result 
places an upper limit on the average deflection suffered by 
a proton in a nuclear collision. It was found that the Intro¬ 
duction of a cos^'fib scattering distribution yielded a satis¬ 
factory result for the gross zenith angle dependence. A 
proton traversing the distance from the top of the atmosphere 
to 3-3 km- will undergo on the average 10 collisions. This 
corresponds to the geometrical cross-section rather than the 
removal path length used above. Its use has been justified 
by Heitler and Janoasy (14). Such protons are extremely rare, 
of course, since they have a path length for removal of only 
twice the distance corresponding to the geometrical cross- 
section. But it makes no difference whether one considers a 
single proton which retains its identity throughout the ten 
collisions or whether one asks for the total angle that the 
final descendant makes with the original primary after ten 
"generations,” the latter situation being far more likely. It 
is permissible, therefore, to consider a single proton. 
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If we assume that each proton suffers exactly ten 
collisions, and that after each collision the new direction 
differs from the old by an angle of exactly + oC, a unique 
distribution after the tenth encounter is determined, which 
depends only on ot . Since oos^'th was introduced artificially 
there is no Justification for asking that this statistical 
distribution coincide with it. One can, however, compare the 
integrals of the two distributions, i.e. ask for the angle 
for which the particle, after ten encounters, suffers a total 
angular deflection of less than, say, twenty degrees, with a 
certain probability computable from the cos^'c b distribution. 
This problem In three dimensions, despite Its similarity to 
the "random walk" problems of mathematical statistics, and 
its simplifying assumptions, is of enormous complexity, an 
unfortunate characteristic of plural scattering problems„ 

The two-dimensional problem, corresponding to the 
projection of the three dimensional scattering distribution 
on a vertical plane, is simpler, and actually corresponds to 
what is observed in this experiment. This two dimensional 
distribution predicts that about 60 % of the particles make 
angles of less than 20 degrees with the primary direction. 

The problem then is to adjust 0 i so that the statistical 
distribution indicates the same thing. The choice of 20 
degrees as the matching point is convenient, for it is close 
to the half value point of the distribtuion. The percentage 
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for a given matching angle is given by another graphical 

integration, hence it is not convenient to match at the exact 

angle within which exactly half the particles are scattered, 

as its determination is unnecessarily laborious. Angles 

giving percentages far removed from 50 # are not expected to 

give the most accurate matching since at small angles the 
6.5 

projected cos k) distribution is too slowly varying, while 
at large angles the statistical distribution is too crude. 

In a single scattering, it is assumed that a particle 
is deflected through an angle of + o< and that the ten 

10 

encounters are statistically independent. Then there are 2 
possible trajectories, corresponding to six distinct total 
angular deflections. These deflections, along with their 
relative probabilities, are enumerated In Table VI. One finds 
that a total scattering angle of about 4oc corresponds to a 
probability of about 60#. Hence from 4oc » 20°, one gets 
oC - 5 °. 

A check of this figure may be had by noticing that a 
proton's behavior is similar to a simple diffusion process. 

The root-mean-square total deviation of a proton undergoing 
N encounters with a scattering angle of OC per encounter is 
given by 

A = Ktf> = fN” oc. 

— o 

With A roughly equal to 20 as before, and with ten encounters. 
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oc = 20 /'Ro ~ 

which compares with the previous result, but the smallness 
of N limits the accuracy of ot . 

Thus a model of the scattering process has been con¬ 
structed which gives a result roughly approximating the 
observed distribution. Present data do not warrant a more 
detailed analysis. 


TABLE VI 

SIMPLE TWO-DIMENSIONAL SCATTERING PROBABILITIES 


e 


Relative Probability Probability That a Particle 
of Scattering Be Scattered Through an 

at Angle $ Angle Less Than 0 


10 

8 

6 

4 

2 

0 


0 

2 

20 

90 

240 

420 

252 


1.0000 

0.9981 

0.9785 

O. 89 O 6 

O .6563 

0.2461 

0.0000 


Total - 1024 - 2 




CHAPTER VI 


COMPARISON WITH THE RESULTS OF OTHER AUTHORS 

The only published data at all concerned with the 
problem treated here is the work of Walker (33) and Branch 
( 36 ). Walker performed a counter-absorber type of experiment 
at 3-26 km. elevation and at sea level in which recorded 
events are those due to charged particles producing penetrat¬ 
ing showers in blocks of lead absorber of varying thicknesses. 
After elaborate precautions to remove the soft and the mesonic 
components, he Is able to call such particles protons, but the 
method of selection only vaguely defines accepted momentum. 
Walker’s requirements for the lowest energy case are that a 
charged particle traverse two trays of counters separated by 
12 inches, enter an 8-inch block of lead, and produce a shower 
there, whose presence shall be indicated by the discharge of 
three or more counters In a tray located below the absorber. 
Knowledge of which counters in the upper two trays were dis¬ 
charged yields the zenith angle dependence to 25 ° in six 
steps. The only limit that can be placed on the extent of 
the accepted momentum band is a lower one, which, because of 
the requirement for shower production, is certainly very much 
higher than the 1 Bev/c maximum cutoff for this experiment. 
This fact is further substantiated by Walker's ratio of the 
rate at elevation to that at sea level. This ratio is about 
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four, which is far removed from the well-established ratio of 
about twenty as obtained by this laboratory, by Wilson (21) 
and by Whittemore and Shutt (3*0* for low (less than 1 Bev/c) 
protons. 

Walker finds that his data at elevation can best be 
correlated with a cos n 0 distribution with n equal to 5-5 + 
0.7, while at sea level the exponent is 3*^ + 1.0. The 
increase in exponent with altitude is not at all reasonable 
in the light of the model presented here. The sea level 
exponent agrees with the figure given in Chapter V, Section 
2, but of course the agreement has little meaning as the data 
do not overlap. Similar work has been done by Brown and McKay 
(37). 

Branch has determined the angular distribution about 

the shower axis of penetrating particles in extensive air 

showers. Such showers, of course, include, in addition to 

nucleons, many pi and mu mesons. A conclusion therefrom, 

relevant to this problem, is nevertheless possible, since 

Branch finds a root mean square angular deviation of emergent 

particles of about 3 degrees. This figure agrees fairly well 
o 

with the figure 5 presented here. These two numbers are to 
be contrasted with the value 15-20 degrees which the early 
theoretical treatments by Messel et al (8-20) gave. Branch's 
result, which this experiment Independently corroborates, has 
caused drastic changes in the nature of the assumed cross- 
section ( 38 , 39 ). 




CHAPTER VII 


SUMMARY 

The proton spectrum for particles in the momentum 
band 0.4-1.0 Bev/c and the meson spectrum for particles in 
the momentum band 0.09-0.34 Bev/c at zenith angles 45° East* 
vertical, and 45° West, have been obtained at 3-3 tan* eleva¬ 
tion and magnetic latitude 48° North. 

The behavior of the zenith angle dependence of mu- 
mesons has been found to be consistent with a transformation 
described by Charbonnier for higher energies, involving only 
the momentum loss by ionization of mesons at minimum ioniza¬ 
tion and the well-defined radioactive properties of this 
particle. 

The zenith angle dependence of protons, along with a 
simplified phenomenological calculation of same, is presented 
as new data on the nucleonic cascade. Theoretical calcula¬ 
tions presently found in the literature are not yet well 
enough founded to allow prediction of these data, nor have 
they been used to yield the differential momentum spectra of 

limited extent as presented here. It is found experimentally 

3.2 

that the gross zenith angle dependence will fit a cos B 
distribution at 45 degrees East or West. The simplified 
calculation assumes that a unidirectional component of the 
primary radiation wills (1) suffer exponential absorption 





along its original path with a path length for removal of 

p 

120 g/cm of air, and (2) acquire an ultimate distribution 
at 3-3 kw- due to scattering fitting a cos n k> law, n being 
an adjustable parameter. 

It is found that a value of n of about 6.5 yields a 
unidirectional distribution which, when integrated over all 
angles, results in a gross zenith angle dependence agreeing 
with cos 0 ' 6 . The cos & distribution with n equal to 6.5 
is quite sharp and leads to the conclusion that the average 
projected scattering angle per collision is only about 5 
degrees. 
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